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Abstract. We study the phase diagram of the two-dimensional fully frustrated XY model 
(FFXY) and of two related models, a lattice discretization of the Landau-Ginzburg- Wilson 
Hamiltonian for the critical modes of the FFXY model, and a coupled Ising-XY model. We 
present Monte Carlo simulations on square lattices L x L, L < 10 3 . We show that the low- 
temperature phase of these models is controlled by the same line of Gaussian fixed points as 
in the standard XY model. We find that, if a model undergoes a unique transition by varying 
temperature, then the transition is of first order. In the opposite case we observe two very 
close transitions: a transition associated with the spin degrees of freedom and, as temperature 
increases, a transition where chiral modes become critical. If they are continuous, they belong 
to the Kosterlitz-Thouless and to the Ising universality class, respectively. Ising and Kosterlitz- 
Thouless behavior is observed only after a preasymptotic regime, which is universal to some 
extent. In the chiral case, the approach is nonmonotonic for most observables, and there is a 
wide region in which finite-size scaling is controlled by an effective exponent v a g w 0.8. This 
explains the result v « 0.8 of many previous studies using smaller lattices. 



1. Fully frustrated systems 

In the last few decades there has been a considerable interest in the consequences of frustration 
on the critical behavior of statistical systems. The simplest example is the antiferromagnetic 
Ising model on a triangular lattice, whose Hamiltonian is 



where J is positive, = ±1, and the sum is extended over all lattice nearest neighbors (ij). 
For T — > the system tends to be antiferromagnetically ordered, i.e. spins on nearest-neighbor 
sites prefer to be oppositely aligned. However, this is not possible everywhere. For instance, 
see Fig. 1, on any lattice triangle one link must be frustrated, i.e. spins on the corresponding 
sites must be parallel so that the local energy assumes its maximum value. The presence of 
frustration has an important consequence. At variance with the ferromagnetic Ising model the 
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Figure 1. A frustrated triangle. If a a = 
1, the local energy associated with links 
AB and AC is minimized by taking gb = 
oc = — 1- The local energy on link BC 
corresponds to a maximum: link BC is 
frustrated. 



antiferromagnetic one is disordered at any temperature: the large entropy forbids an ordering 
transition [1, 2]. 

The Ising model can be generalized by considering the N- vector model on a triangular lattice. 
In this case one considers unit A^-component vectors Sj and the Hamiltonian 



Also this model is frustrated: There is no configuration in which all neighboring spins are 
antiparallel. However, at variance with the Ising case, here the entropy vanishes at zero 
temperature. Indeed, once rotational invariance has been broken by fixing the direction of one 
spin, there is a finite number of configurations that are global minima of the Hamiltonian [6]. 
For N = 2, the only case we will consider, if s*j = (cos#j,sin0j), one must have \9i — 0j\ = 2ir/?> 
or 47r/3 when i and j are nearest-neighbor sites. It is easy to verify that the degeneracy of the 
ground state is 7Li % 0(2), where 0(2) is the invariance rotation group. The group Z2 is due to 
the possibility of two (chirally) different configurations. As shown in Fig. 2, the ground state 
is uniquely determined once one breaks rotational invariance (by setting, for instance, 9 a = 0) 
and chooses the chirality of triangle ABC (by setting 9b = 120° or 240°). An observable that 
distinguishes between the two ground states is the chirality. Given a lattice triangle, see Fig. 2, 
we can consider [6] 



which assumes the values ±1 on any lattice triangle in a ground-state configuration. A good order 
parameter is obtained as follows. We assign s n = ±1 to each lattice triangle so that s n = —s m if 
triangles n and m share a lattice link. The order parameter, the chiral magnetization, is simply 




(2) 



(ij) 



C n = —{ S m(9 A - Ob) + sin(0 B - 0c) + sin(0 c - 0a)], 



(3) 




(4) 



n 



where the sum is extended over all lattice triangles. 



4 It is interesting to note that this is not true for the spin-5* antiferromagnetic Ising model if S is large enough; 
see [3, 4, 5] and references therein. 
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Figure 2. Two inequivalent ground states related by a chiral transformation. They are obtained 
as follows: one first fixes 6a = 0°, breaking rotational invariance. Then, there are two possible 
choices: on the left we choose 9b = 240°, 6q = 120°; on the right we make the opposite choice. 
All other lattice spins are univocally defined. 



Figure 3. The couplings jij in Hamilto- 
nian (5): jij = 1 on thin lines, jij = —a on 
thick lines. 



It is possible to define a frustrated model also on the square lattice. The relevant model is not 
the antiferromagnetic Ising or iV-vector model, since no frustration occurs on the square lattice 
or, in general, on any bipartite lattice. To obtain a frustrated model we consider a Hamiltonian 
of the form [6] 

T^ffxy = - J y^Jij Si • Sj, (5) 

(ij) 

where the two-component spins Si satisfy Si-Si = 1, jij = 1 along all horizontal lines, while along 
vertical lines ferromagnetic jij = 1 and antiferromagnetic jij = —a (a > 0) couplings alternate, 
see Fig. 3. This model is frustrated for any positive a. Maximal frustration is obtained by taking 
a = 1; for this reason this particular model, the only one we shall consider in the following, is 
called fully frustrated XY (FFXY) model. The square-lattice FFXY model admits two chirally 
different ground states, see Fig. 4, and thus it has the same ground-state degeneracy of the 
antiferromagnetic model on the triangular lattice. These two models are particular examples of 
a general class of systems that all have a Z2<8>0(2) ground-state degeneracy. We will collectively 
call them FFXY systems. 

Even though the symmetry of the FFXY systems (2) and (5) is the same as that of the 
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Figure 4. Ground states of the square-lattice FFXY model: in this case nearest-neighbor spins 
must satisfy 9i — 9j = 45° or 315° if they are connected by a ferromagnetic link, and 9i — 6j = 135° 
or 225° if they are connected by an antiferromagnetic link. Once we fix Op = 0°, there are two 
(chirally) inequivalent possibilities: (left) 6a = 45° or (right) 6b = 45°. All other spins are 
fixed. Note that the ground-state configurations are invariant under translations of two lattice 
spacings. 



ferromagnetic XY model, we expect the critical behavior to be different. Indeed, the universality 
class is not only determined by the symmetry of the order parameter but also by the symmetry 
breaking pattern that is different in the two cases. 

In order to determine the critical behavior we can perform a direct numerical study of the 
model. There is however another possibility, which is the basis of the field-theoretical approach 
to critical phenomena. In this case one first identifies the critical modes of the microscopic 
Hamiltonian and then writes down an effective coarse-grained (continuum) Hamiltonian for 
them. The model one obtains is no longer frustrated; still, it is expected to have the same 
critical behavior as the original one. In order to derive the effective theory, let us consider 
the antiferromagnetic model on a triangular lattice. As is evident from Fig. 2, in the ground 
state spins rotate by 120° when moving in the x direction from one site to its neighbor. Thus, 
critical modes are associated with fluctuations close to the complex Fourier component s(Q) 
with Q = [27r/(3a), 0], where a is the lattice spacing. These fluctuations are parametrized by a 
complex vector <5 = (f)\ + i(f>2 • Note that the appearance of two real two-component fields is at 
variance with the ferromagnetic case. Indeed, in that case, the relevant modes are associated 
with the zero-momentum component s(q = 0). As a consequence of the reality condition 
s(q) = s*(—q), fluctuations are real and are parametrized by a single real two-component field. 
A standard calculation [7, 8, 9] gives the effective Hamiltonian for the fields 4> a : 



^ a=l,2 ' ' a=l,2 



(6) 



A similar argument applies to the square-lattice FFXY model and gives the same effective 
Hamiltonian (6). 

Hamiltonian (6) has a larger symmetry than the original one. Indeed, the symmetry group 
is [0(2) © 0(2)] <g> Z2: the 0(2) groups are related to independent rotations of the two fields, 
while the %2 group corresponds to the field-interchange symmetry. Nonetheless — and this is 



the only property that matters — Hamiltonian (6) for v o > and FFXY models have the same 
symmetry-breaking pattern. Indeed, for vq > 0, the ground state corresponds to (j>\ = and 
4>2 7^ or the opposite and thus it has the same ground-state degeneracy: the ground state 
breaks one of the two 0(2) groups and the Z2 field-interchange symmetry. 
A lattice discretization of (6) is [10] 

= - j e $ a>i • ki + E Ki + u «i - x ) 2 ] + 2{u + D ) E A> ( ? ) 

(ij'),o a,i i 

where J > (the model is ferromagnetic), a = 1,2, <^> aj j is a real two-component variable, the 
first sum goes over all nearest-neighbor pairs, and = cj) a ■ (f) a . The correct symmetry-breaking 
pattern is obtained for D > 0, which corresponds to vq > in (6). Moreover, stability requires 
U > 0. For U — > 00, 7^0 becomes simpler and we obtain 

H = - J E 4i • + 2D E ( 8 ) 

<ij"),o i 

where the fields satisfy the constraint (fif i + <^>| i = 1. This is the 4- vector model with a spin-4 
perturbation that breaks the 0(4) symmetry to [0(2) © 0(2)] <g> Z2. If we additionally take 
D — > 00 we must have 4>\ ^ j = 0. In this case, we can parametrize 

<j>i,i = |(1 + cri)sj, 02,i = 5(1 - (9) 

where cjj is an Ising spin and Sj is a unit two-component vector. The Hamiltonian reduces to 

W = -|X)(l + ^)ai-sj-. (10) 

<«> 

Hamiltonian (10) has the same invariance as (7) although, in terms of the new fields, the 
0(2) © 0(2) symmetry is nonlinear ly realized: 



Si = [1(1+^1 + ^(1-^)^2]^ (11) 



< = <*» (12) 

where R\ and R2 are 0(2) rotation matrices. It is possible to add terms of the form a^Gj without 
breaking the symmetry of the Hamiltonian. We can thus consider the more general Hamiltonian 

[11] 

-(1 + GiCFj) Si ■ Sj + CUiUj 



Ti-isXY = - E 

(ij) 



(13) 



We will call this model the Ising-XY (IsXY) model. For J > and C + J/2 > it has the 
same symmetry-breaking pattern as FFXY systems. Thus, it represents another example of this 
class of models. We should mention that there are many other systems that share the same 
symmetry-breaking pattern: for an extensive list of references see [12]. 



2. Results 

Two-dimensional FFXY systems (but note that these systems have also been studied, both 
theoretically and experimentally, in three dimensions [13, 14, 15, 16]) have been extensively 
studied in the last thirty years, after the appearance of the seminal papers by Villain [6]. For an 



extensive list of references, see [12]. In spite of that, their critical behavior is object of debate 
still today. Two scenarios have been proposed for the critical behavior of models (2) and (5). 

A first possibility is that these models have two continuous transitions. As temperature 
decreases, there is first a transition associated with the chiral degrees of freedom: at the 
transition there is no magnetic ordering but only chiral order. As temperature further decreases, 
there is an intermediate phase in which spins are disordered while chiral variables are magnetized. 
Then, a second transition occurs, followed by a low-temperature (LT) phase in which spin-spin 
correlations decay algebraically. In this scenario chiral and spin modes do not interact at the 
transitions and thus, if the transitions are continuous, one expects a chiral Ising transition (the 
order parameter is a scalar) and a spin Kosterlitz-Thouless (KT) transition. This scenario has 
not been confirmed numerically so far. The computed exponents at the chiral transition do not 
agree with the Ising ones. For instance, one finds v rs 0.8 instead of the Ising value v = 1. 
Moreover, it is not clear how much one can believe in the presence of two transitions that are 
very close to each other: numerically (T sp i n — T c hi ra i) — 10~ 2 J. 

The inconsistencies of the two-transition scenario apparently favor the presence of a single 
critical point. In this scenario, the observed difference between the critical temperatures is 
interpreted as a correction-to-scaling effect. Since chiral and spin modes become critical at the 
same temperature, it is possible that the transition belongs to a new universality class with a 
new set of critical exponents. In this scenario, the result v « 0.8 would be fully acceptable. Of 
course, it is also possible to interpret the results for v in terms of crossover effects. Since the 
two transitions are very close, scaling corrections may be large so that the asymptotic behavior 
can be observed only on very large lattices [17]. 

We have considered again the issue [10, 12], performing extensive simulations of the FFXY 
model (5), of the 4 model (7) with U = 1 and of the IsXY model (13). In all cases we have 
considered the square lattice. We have used a mixture of Metropolis and overrelaxation updates 
as well as cluster updates in the LT phase, essentially following [18] (see [12] for a detailed 
discussion). We have studied the critical behavior on lattices of size L x L, in some cases up to 
L ~ 10 3 : for the FFXY model, the 4 model with D = 1/2, and the IsXY model with C = 0, 
the largest lattice we have used at the chiral transition corresponds to L = 1000, L = 1200, 
L = 360 respectively. In the LT phase the Monte Carlo algorithm is much more efficient and 
we have been able to simulate even larger sizes: for the <fi A model with D = 1/2 we performed 
simulations for L = 2048. 

The analysis of the Monte Carlo results for the square-lattice FFXY model definitely shows 
that this model undergoes two transitions: the chiral one belongs to the Ising universality class, 
while the spin one is compatible with a KT behavior. In the LT phase the critical behavior of 
the spin modes is controlled by the same line of Gaussian fixed points as in the standard XY 
model. The discrepancies from Ising behavior at the chiral transition that have been observed in 
previous studies are simply crossover effects. They are due to the presence of a large, albeit finite 

(c) 

spin correlation length Q at the chiral transition. In finite-size scaling studies the asymptotic 
behavior can only be observed if L S> ffi . Since ffi is quite large, £^ = 118(1), the asymptotic 
behavior can only be observed in simulations with L ~ 500-1000, i.e. for values of L that 
are much larger than those that could be used in simulations until a few years ago. In the 
other models we have studied, the determination of the asymptotic behavior may be even more 
difficult. For instance, in the (j) 4 model with D = 1/2, £g C ^ ~ 380. Thus, even with simulations 
with L = 1200, we have not been able to observe the Ising behavior but only the beginning of 
the crossover towards the asymptotic behavior. 

In the IsXY and in the <fi 4 model we have observed two transitions in a large parameter 
region. However, for D or — C large, we have also observed a unique first-order transition which 
separates the LT phase with chiral order and spin quasi-long-range order from the disordered 
phase, see Fig. 5. Thus, our results confirm the two-transition scenario for generic FFXY systems 




Figure 5. Sketch of the phase diagram of the (j) 4 model (7) for U = 1 and D > (left) and 
of the IsXY model (13) (right). The continuous, dashed, and thick continuous lines represent 
Ising, KT, and first-order transition lines. The distance between the ferromagnetic Ising and 
KT lines is amplified; otherwise, the two transitions cannot be distinguished on the scale of the 
figure. The phase diagram within the circled region is unknown. In the IsXY case there is also 
an antiferromagnetic Ising transition (af) starting at C = —C\ s = — 1(1 + \/2), J = 0. 

in the sense that we have found no evidence of a unique continuous transition where chiral and 
spin modes become both critical. A single transition occurs only if it is of first order. 

Beside these results that confirm the two-transition scenario, we have also observed 
an unexpected universal crossover behavior. We find that renormalization-group invariant 
quantities (e.g., critical exponents, Binder parameters, . . .) computed in the different models 
scale at the chiral and spin transitions respectively as 

K = f£\L/i), K = f£\L/l), (14) 

where I is a model-dependent scaling factor that is identical at the two transitions. At the 
chiral transition — the only case in which we have done a systematic investigation by varying C 
and D — corrections appear to increase as D or — C increases (in practice, significant deviations 
are observed for D > 4 and C < — 2). Of course, it is of interest to have a renormalization- 
group explanation of this apparent universality. Eq. (14) can be explained by the presence of 
a multicritical point [19, 20, 21] or a line of multicritical points of the same type) where chiral 
and spin modes become both critical. Indeed, close to a multicritical point we expect that any 
RG-invariant quantity behaves as 

n = f n (L/^L/i ch ), (15) 

where £ s and £ c h are the infinite-volume correlation lengths for spin and chiral variables. Our 
analysis at the Ising and KT transitions corresponds to fixing L/£ c h = and L/^ s = 0, i.e. 
provides the scaling function along two particular lines. If the interpretation in terms of a 
multicritical point is correct, the functions fii(x) provide informations on the behavior at the 
multicritical point. Indeed, while Ising or KT behavior is observed for x — > oo, in the opposite 
limit i->0we obtain the value of the RG-invariant quantity 1Z at the multicritical point. 

The nature of the multicritical point is unclear. One possibility is the 0(4) multicritical 
point that is present in the </> 4 theory for D = 0. Another possibility is the multicritical point 
that appears in frustrated XY systems with modulated couplings (for instance in model (5) for 
a / 1) [22, 23, 24] or in generalizations of the IsXY model, in which an additional spin-spin 
coupling is added, breaking the 0(2) © 0(2) symmetry. 

Finally, we wish to compare with field-theory (FT) approaches. Perturbative analyses [25] 
of model (6) indicate the existence of a new universality class associated with the symmetry- 
breaking pattern [0(2) © 0(2)] © Z2 — > 0(2). Even though we have found no evidence for it, 
our results do not necessarily contradict those of Ref. [25]. It is possible that the models we 



have considered are outside the attraction domain of the FT fixed point. If this is the case, field 
theory provides another candidate for the multicritical point. The models we consider could be 

(c) 

outside, but close to the attraction domain of the fixed point — this is not unplausible since Q 
is large — so that the crossover behavior is controlled by the FT fixed point. We wish also to 
make a remark on the validity of (6). In the derivation of Hamiltonian (6) by using the standard 
Hubbard-Stratonovitch transformation, terms with more than four fields are neglected [7, 8]. 
In particular, terms of the form (0i • 4>2) n appear at sixth order (n = 3) (resp. eighth order) 
in the case of the triangular-lattice (resp. square-lattice) FFXY model [8]. These terms have 
only Z2 © 0(2) symmetry and thus, under renormalization-group transformations, are bound to 
generate a term of the form (0i • (P2) 2 , or even a quadratic term of the form (</>i • ^2). We obtain 
therefore the multicritical Hamiltonian 



Wlgw,2 = Hlgw + J d d x 



\r 2 {$i ■ $2) + \zSi ■ 4>if + \z\$\ ■ fa){<t>l + <&) 



(16) 



as it has been postulated for modulated systems. If this interpretation is correct, the FT fixed 
point may only be relevant at the multicritical point, provided it is stable under the quartic 
perturbations that break the 0(2) © 0(2) symmetry. This holds in three dimensions [26], but 
nothing is known in the two-dimensional case. It should be remarked that these considerations 
are only relevant for the FFXY model. The IsXY and 4> 4 theories are 0(2) ©0(2) invariant and 
thus the correct FT Hamiltonian is clearly (6) and not (16). 



3. Chiral transition 

In order to determine the nature of the chiral transition, we have studied the behavior of several 
quantities at fixed R c = £ C /L where £ c is the chiral correlation length (see [12] for a precise 
definition in the different models). We use the method proposed in [27] and further discussed in 
[28]. We fix R c equal to i? Is , where Ri s = 0.9050488292(4) is the universal value of £/L at the 
critical point in the 2-d Ising universality class [29]. We stress that this choice does not bias our 
analysis in favor of the Ising nature of the chiral transition. For any chosen value (as long as it is 
positive) and whatever the universality class of the chiral transition is (it may also coincide with 
the spin transition), we are studying the model for L-dependent temperatures T e g (L) such that 
T e ff(L) — > T c h for L — > 00. Note, however, that quantities like the Binder parameter depend 
on the chosen value for R c . Indeed, in the finite-size scaling limit, R c = /i?[L 1 / jy (T — T c h)]. 
Therefore, fixing R c is equivalent to fixing X = L 1 ^ (T c r(L) — T c h). Since the Binder parameter 
satisfies an analogous relation B c = /£[L 1 /"(T — T c h)], at fixed R c B c converges to fs(X). By 
fixing R c to the critical Ising value, we will be able to perform an additional consistency check. 
If the chiral transition belongs to the Ising universality class, then X = (apart from scaling 
corrections) and we should find that any RG-invariant quantity converges to its critical-point 
value in the Ising model. 

We first verify that £ s converges to a constant as L — > 00. In Fig. 6 we show the numerical 
results. The correlation length is clearly finite in the FFXY model and in the IsXY model with 
C = 0. In the 4> 4 model with D = 1/2 we do not yet observe that £ s is finite, although it is already 
clear that £ s does not increase linearly with L, as it would be the case if the spin correlation 
length were infinite for L = 00. Then, we verified that the transition, if continuous, belongs to 
the Ising universality class. The best evidence is provided by the Binder chiral parameter. If 
the transition belongs to the Ising universality class we should find [29, 30, 12] 

B c = B Is + bL- 7 /\ (17) 

where B\ s = 1.167823(5) [29]. The results reported in Fig. 7 are fully consistent with Ising 

(c) 

behavior for L 3> Q ■ Not only do we observe the Ising asymptotic value, but also the rate of 
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Figure 6. Spin correlation length £ s at fixed R c = R\ s (chiral transition). For L — > oo, we have 
^ = 118(1) in the FFXY model, £ s = 52.7(4) in the IsXY with C = 0, and £ s « 380 in the c/> 4 
model with U = 1 and D = 1/2. 



0.04 r 

0.03- 

AB 

c 

0.02 
0.01- 



1000 



.i - 



0.0001 



I 

* 



• FFXY 
■ lsXYC=0 
▼ 4>, D=l/2 



-7/4 



0.0002 



Figure 7. Chiral Binder parameter B c at fixed i? c = R\ s (chiral transition). Plot of 
AB C = B c - Bi s at fixed R c = Ri s vs L" 7 / 4 , for the FFXY model, the c/> 4 model at D = 1/2, 
and the IsXY model at C = 0. Bi s = 1.167923(5) is the value of the Binder parameter at the 
critical point in the Ising model [29]. 



convergence is well verified. Also the critical exponents v and rj converge to the Ising values for 
L S> £^ (results for v will be shown below). 

4. Low-temperature phase and spin transition 

Since the spin correlation length is finite at the chiral transition there must be a paramagnetic 
phase with chiral order. Such a phase ends at a second transition which is followed by the LT 
phase in which chiral order and spin quasi-long-range order coexist. 

We first study the nature of the LT phase and verify the breaking of the Z2 invariance. 
Direct evidence is provided by the chiral Binder parameter B c . If chiral modes are magnetized, 




Figure 8. Estimates of R s = £, S /L y s rj in the LT phase. The continuous line is the prediction 
obtained by assuming that in the LT phase criticality is controlled by the same line of Gaussian 
fixed points as in the XY model. 



B c — ► 1 + 0(L~ 2 ) for L — > oo. This behavior is very well verified in all models. In the 
model, the Z2 group corresponds to the field-interchange symmetry. Thus, the Z2 symmetry 
breaking implies that only one of the two fields 4>i an d 4>2 is critical in the LT phase. Our 
numerical results fully confirm this expectation. One can distinguish the fields according to 
the value of Q a = X^^a- The field with the largest value of Q is critical (for instance, the 
corresponding susceptibility and correlation length diverge as L — > 00), while the other one is 
not (the susceptibility and the correlation length have a finite limit as L — > 00). 

Once we have checked that chiral modes are magnetized (this is of course obvious because of 
the presence of the chiral transition), we study the behavior of the spin variables and verify that 
the large-L behavior is controlled by the same line of Gaussian fixed points as in the standard 
XY model. In the XY model one can derive universal relations among renormalization-group 
invariant quantities that are valid in the whole LT phase, up to the KT transition [31]. Indeed, 
below the KT transition the spin-wave approximation is asymptotically exact as L — > 00 and 
allows an analytic determination of any quantity in terms of the spin- wave parameter. Such a 
parameter is not universal and can be eliminated by expressing a renormalization-group invariant 
quantity in terms of another. For instance, one can express the helicity modulus on a square 
lattice L x L with periodic boundary conditions in terms of the exponent rj computed from the 
size dependence of the magnetic susceptibility (\ ~ L 2_,? ): 

T = _L_ Y.n=-oc n2 eM-™ 2 h) 

exp(-7rn 2 /?7) ' 

where < 77 < 1/4. Analogously, one can express £ S /L in terms of 7]. In Fig. 8 we compare 
the spin-wave prediction with numerical results. The agreement is quite good, confirming that 
FFXY systems and the standard XY model have the same LT phase, as far as the spin degrees 
of freedom are concerned. 

The above-reported results for the LT phase make it plausible that the spin transition belongs 
to the KT univerality class. Another check is provided by our numerical results for £ S /L and 
T. In the XY model, at the KT transition [31] £ S /L « 0.750691 + 0.212430/ ln(L/C x ) and 
T w 0.636508 + 0.318899/ ln(L/C2). In all models we have studied these two quantities assume 
the XY values approximately at the same temperature — which we identify with the spin critical 
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Figure 9. Ratio R s = £ S /L (£ s is the spin correlation length) at the chiral transition vs L r = L/l 

for the FFXY model and 4> 4 (left) and IsXY (right) models. We set I = for the FFXY model; 
the values of I for the other models are obtained by requiring that all data fall on a single curve. 
Note the logarithmic scale on both axes. For L r — * oo, R s converges to 0. 



temperature — thereby confirming the KT nature of the transition. 
5. Crossover behavior 

As we already mentioned, our data show the scaling behaviors (14). In this section we shall give 
a few details. 

In order to verify the scaling behavior (14) at the chiral transition 5 we have first considered 
the data for R s = £ S /L (£ s is the spin correlation length) at the chiral transition and we have 
investigated whether they fall on a single curve by using a rescaled variable L r = L/l, where / is 
a rescaling factor that depends on the model. The result are reported in Fig. 9. The data fall on 
a single curve with remarkable precision, i.e. £ S /L = f s (L/l) where f s (x) is model independent. 
Note that the rescaling factors change significantly from one model to another: for instance, 
I/Iffxy = 0.75 (resp. 7.0) in the (p 4 model with D = 4 (resp. D = 1/5) and Z/Z F fxy = 0.031 
(resp. 12) in the IsXY model with C = 0.3 (resp. C = —2). It is easy to realize that I should be 

(c) 

proportional to Q , the infinite- volume spin correlation length at the chiral transition. Indeed, 
since £ s — > £s as L — ► oo, we have f s (x) ~ a/x for x — ► oo, where a is model independent (f s {x) 
is model independent). Moreover, = l/a. Therefore, if we fix I = Q c ' in one model, then 

(c) 

the same holds in all different models. In the figures we have chosen Zffxy = 118 ~ Q and 
thus the plots we present are indeed in terms of L/Q c \ even in those cases in which we have 
not been able to determine directly the spin correlation length for L — > oo. 

In order to verify the universality of the scaling Ansatz (14), we have considered the spin 
Binder parameters B S( p and B s (we use two different inequivalent definitions, see [12]). In Fig. 10 
we plot B s and B s ^ versus L/l. We use the rescaling factors that have been determined in the 
analysis of R s . The agreement is quite good. Deviations appear as D or — C increases. In 
particular, the data for D = 4 (cp 4 model) and for C = —3 (IsXY model) are outside the curve. 
They would fall on the same curve as the others, only if the rescaling factor is changed by a 
factor of 2 and 5 respectively in the two cases. In Fig. 11 we plot the results for the helicity 

Note that our results have been obtained at fixed R c and not at the chiral critical point. However, since we are 
dealing with an Ising transition and R c has been fixed to the critical-point Ising value, this is irrelevant in the 
scaling limit. Had we fixed R c to a different value, we would have obtained quantitatively different scaling curves, 
corresponding to the limit £ c h,L — > oo at fixed L/£ch 7^ in Eq. (15) (£ c h is the infinite- volume chiral correlation 
length) . 
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Figure 10. Spin Binder parameters B s and B s ^ at the chiral transition vs L r = L/l. Results for 



the FFXY, </> 4 , and IsXY models. For L r - 
The rescalings / are the same as in Fig. 9. 



oo, B s and converge to 2 and 3/2 respectively. 
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Figure 11. Helicity modulus T (left) and chiral Binder parameter B c (right) at the chiral 
transition vs L r = L/l. We report AB C = B c — Bi s , where B\ s is the value of the Binder 
parameter at the critical point in the Ising model. The rescalings I are the same as in Fig. 9. 



modulus: again all data fall on a single curve quite precisely. It is interesting to note that, for 
0.02 < L r < 0.5 R s and T show an approximate power-law behavior: they behave as L~ e with 
e ~ 0.1 (R s ) and e ~ 0.33 (T). If this behavior holds also for smaller values of L r , we have 
T, R s — » oo for L r — > 0. 

Next we consider the chiral variables. They also show the universal behavior (14). In Fig. 11 
we report the Binder chiral parameter, with the same scaling factors I as before. The data fall 
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Figure 12. Effective exponent 1/V e ff computed by using R c = £ C /L, the Binder chiral parameter 
B c , and R s = £ S /L (£ c and £ s are respectively the chiral and spin correlation lengths) vs L r = L/l 
at the chiral transition. For L r — ► oo, l/^ e ff converges to l/v\ s = 1 for R c and B c , and —1 for 
R s . The rescalings Z are the same as in Fig. 9. 



again on a single curve although the curve does not change significantly as L/l varies, at variance 
with the spin variables. 

Finally, we show the effective exponent v e R that can be obtained from the derivatives of 
R c , B c , and R s . The exponent obtained from the chiral variables should converge to the Ising 

value v = 1 as L — > oo, and indeed it does. The approach is however nonmonotonic, u e g being 

(c) 

first smaller than 1, then larger. It is interesting to note that for L r < 1, i.e. L < Q , v e g is 
approximately constant and equal to 0.8. This behavior explains previous results. Indeed, if 
one performs simulations only for values of L such that L r < 1 (this was the case in previous 
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Figure 13. Ratio R s = £ S /L (left) and helicity modulus T (right) at the spin transition vs 
L r = L/l. The rescalings I are the same as in Fig. 9. 



simulations of the FFXY model since 

d c) « io 2 ) one would estimate v = 0.8. Note also that, 

tc) 

for L < <Q , the effective exponent f e fT obtained from the spin variable R s is approximately 
constant and close to the value obtained by using chiral variables. In this range of values of L 
chiral and spin variables appear as if they are both critical. 

The same analysis can be repeated at the spin transition. In Fig. 13 we report £ S /L and 
T versus L/l where the rescaling factors I are those determined at the chiral transition. The 
agreement is again quite good. The existence of scaling at the two transitions with the same 
rescaling factors is another piece of evidence in favor of the multicritical origin of the universality 
we observe. 

The crossover curves we have computed can give us some hints on the nature of the 
multicritical point, if it really exists. Indeed, the behavior at the multicritical point is simply 
obtained by considering the limit L r — > 0. First, let us notice that the data at the chiral 
transition apparently exclude the possibility of a decoupled multicritical point in which spin and 
chiral modes have XY and Ising behavior. Indeed, the helicity modulus is very much different 
from the KT value, T kt = 0.63650817819. . . [31] and B c is apparently smaller than the Ising 
value for L r — * 0. 0(4) behavior is possible, since, as we already discussed, our data are 
compatible with T,R S — > oo for L r — > 0. Also the data for the Binder parameters B s and 
are compatible with the 0(4) value, B s = B st f, = 1. For L r — > 0, the crossover curves at the spin 
transition should converge to the same values as those at the chiral critical point. This is not 
evident from the results plotted in Fig. 13. This can be easily explained. At the spin transition 

(s) 

the natural scale is the infinite- volume chiral correlation length at the transition £c • We expect 
XY behavior for L/£c 3> 1 and multicritical behavior in the opposite case. Numerically, we 

find £s /£c ~ 15, so that L/£c = 1 corresponds to L r = L/^ « 0.07. As it can be seen from 
the figure, none of our data satisfies the condition L r <C 0.07, so that we are unable to observe 
multicritical behavior at the spin transition. 
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